Optical precursors in transparent media 
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We theoretically study the linear propagation of a stepwise pulse through a dilute dispersive 
medium when the frequency of the optical carrier coincides with the center of a natural or elec- 
tromagnetically induced transparency window of the medium (slow-light systems). We obtain fully 
analytical expressions of the entirety of the step response and show that, for parameters repre- 
sentative of real experiments, Sommerfeld-Brillouin precursors, main field and second precursors 
("postcursors") can be distinctly observed, all with amplitudes comparable to that of the incident 
step. This behavior strongly contrasts with that of the systems generally considered up to now. 

PACS numbers: 42.25.Hz, 42.25.Kb, 42.25.Lc 



As far back as 1914, Sommerfeld and Brillouin theoret- 
ically studied the propagation of a stepwise pulse through 
a linear dispersive medium [l|, 3, They showed in 
particular [2| that the arrival of the main signal is pre- 
ceded by that of two successive transients they named 
forerunners. The first one (now called the Sommerfeld 
precursor) arrives with the velocity c of light in vacuum. 
Its instantaneous frequency, initially higher than the fre- 
quency uJc of the optical carrier, decreases as a function 
of time whereas that of the second one (the Brillouin pre- 
cursor), initially lower than ojc, evolves in the opposite 
direction. Sommerfeld and Brillouin considered a single- 
resonance Lorentz medium and made their calculation 
by using the newly developed saddle-point method of in- 
tegration. Revisited by various methods, this problem 
has become a canonical problem in electromagnetics and 
optics [1, H, Different models of medium have obvi- 
ously been considered and the theoretical literature on 
precursors is very abundant. See [3| for a recent review. 

As intuitively expected, the precursors will be observed 
only if the rise-time of the incident step is short com- 
pared to the response time of the medium Q. Most of 
the theoretical papers consider dense media with very 
short response time (< 1 fs) and the fulfillment of the 
previous condition raises serious experimental difficulties. 
This explains the dramatic dearth of papers reporting di- 
rect demonstrations of precursors. A first experiment was 
achieved in the microwave region with waveguides whose 
dispersion mimics that of the Lorentz medium ^j. In the 
optical domain, Aavikssoo et al. studied the propaga- 
tion of sing le-ended exponential pulses through a GaAs 
crystal [10||. Associated with an exciton line, the precur- 
sors then appear_as a spike superimposed on the main 
pulse (see also [11]). A discussion on the observability of 
optical precursors in dense media can be found in |12]. 

Much more favorable time scales are obtained by ex- 
ploiting the narrowness of atomic or molecular lines in 
vapors or gases. The switching times of the incident field 
may then be very long compared to the optical period 



without washing out the transients. In such conditions, 
the slowly varying envelope approximation (SVEA) is ab- 
solutely justified. The medium is fully characterized by 
its system function 77(0) connecting the Fourier trans- 
forms of the envelopes of the transmitted and incident 
fields . r2 designates the deviation of the current op- 
tical frequency lo from the carrier frequency ujc and the 
envelope of the optical step response reads 



a{t) 



(1) 



where the contour F is a straight line parallel to the real 
axis passing under the pole at = 0. Eq.|(l]) can always 
be numerically solved by means of fast Fourier transform 
(FFT) but, generally, has no analytical solution. For- 
tunately enough, such a solution exists in the reference 
case of a medium with a single Lorentzian absorption- 
line (see, e.g., [3|)- On resonance and for large optical 
thickness, a{t) takes the simple form 



a(t > 0) = e"^* Jo (V2q^) 



(2) 



where L is the medium thickness, t (as in all the fol- 
lowing) is a local time (real time minus L/c), a is the 
resonant absorption-coefficient for the intensity (a/2 for 
the amplitude) and 7 is the half width at half maximum 
of the line. For t > ti = jijTj^ the asymptotic form of Jq 
may be used and a{t) approximately reads 



a{t > ti) 



-^t cos (72^1^ -7r/4) 



Experimentally evidenced in [13|, the transient given by 
Eq.([2|) and Eq.|[3]) may be formally analyzed in terms of 
Sommerfeld and Brillouin precursors, which are tempo- 
rally superimposed in dilute media [la, [l3|- However we 
remark that these "precursors" precede nothing since the 
medium is then opaque for the "main field". In order to 
obtain true precursors we examine in this letter the much 
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richer case where the medium is (nearly) transparent at 
ujc ■ Our main purpose is to establish aproximate analyt- 
ical expressions of the step response of such media, FFT 
being used to check the validity of the approximations. 

We consider first a medium with a natural trans- 
parency window between two identical absorption-lines of 
intensity optical thickness aL/2 ^ 1 located at ± A- 
Such a medium has proved to be a very efficient slow- 
lightsystem (isl . [l9l . Eol . [2l| . Its system function reads 

mm 



Eq.JT])) reads 



H{n) = exp 



aL"f 



1 



1 



^ + i(ri + A) 7 + z (17 - A) 



(4) 

A good transparency at = is achieved if 7 ^ 
A and aLj'^/A'^ < 1. The group delay then reads 
Tg = aL7/2A2 [2^1 and H{0) = exp (-aL7V2A2) = 
exp(— 7T£,). Figd] shows the step response obtained for 
parameters representative of the slow-light experiments 
achieved on a cesium vapor in the near infrared The 




1 2 
Time in group-delay units 

Figure 1: Step-response a(t) of a medium with a natural 
transparency window. The analytical (•) and numerical (full 
line) forms are respectively obtained by asymptotic calcula- 
tions (see text) and by the means of a FFT involving 2^^ 
points with a time resolution of 1.2 • lO'^Tg. The step of 
amplitude H(0) retarded by Tg is given for reference (dotted 
line). Inset : enlargement of the precursors. The parameters 
are A = 28.9 ns"\ 7 = 0.0164 ns~^ and aL = 2 • 10^ leading 
to Tg = 1.96ns, H{0) = exp(-7rg) = 0.968 and b = 5.22ns-\ 



analytical form is obtained by taking advantage of the 
large value of aL. We note first that, in its very far 
wings, H{Q) equals the system function of a medium 
with a single line of intensity optical thickness aL and, 
as expected, the short time behavior of a{t) is well de- 
scribed by Eq.©. For t > ti ^ = Ia^' can 
be entirely calculated by the saddle point method [23 . [25| . 
The significant contributions to a[t) originates in the rel- 
evant saddle points and, eventually, in the pole at = 0. 
Introducing the phase function ^'(17) = i^t + In [H (fi)]. 



a{t) 



exp[^' {n)]dn/2iTTn 



(5) 



The integral is calculated by deforming F in a contour F' 
traveling along fines of steepest descent of the function 
*(17) from the saddle points where ^''(f^) = . The 
contribution of a non degenerate saddle point at Vis to 
the integral reads 

as{t) = (ins^2^\W'{ns)\) exp [* {^s) + (6) 

where 9s is the angle of the direction of steepest descent 
with the real axis. Note that the instantaneous frequency 
of as(t), defined as d(lTa"^) jdt , equals Re (fig). 

In the present problem, the equation ^''(ri) — giv- 
ing the saddle points can be reduced to a biquadratic 
equation with exact analytic solutions. The latter can be 
regrouped in two pairs i7^(i) = i7± r2„ (i) with n — 1,2 
and 



l-(-l)" Jl + 8t/ 



2t 



(7) 



At every time, Vli (t) is real and very large compared to 7, 
decreasing from Ayjrg/t for i ^ to A for i ^ 00. The 
corresponding saddle points are always non-degenerate 
and their contribution ai [t) = at {t) + [t) to a{t) is 
easily derived from Eq.® with 0^ = ±7r/4. It reads 



ai{t) 



TT 



7* 



cos {^it + A'^TgVtil {^1 - A2) - 7r/4} 



r^iA 



(r!i + A)"^ + (r!i-A)" 



(8) 

As expected, ai(t) tends to a{t) given by Eq.(l3]) when 
t <^Tg. More generafiy, ^2 is purely imaginary for t < Tg 
and the contribution of the corresponding saddle points is 
negligible, except in the vicinity of Tg. So, in a wide time- 
domain, ai{t) is actually the only significant contribution 
to a{t). The corresponding optical field reads Ei{t) — 
Et{t) + E^{t) where Ef{t) = Re [a^{t) exp [ii^Jct)] have 
instantaneous frequencies uof' [t) = uic ± fii {t). Due 
to the time dependence of these frequencies, Ei{t) and 
Ei{t) may be identified respectively to the Sommer- 
feld precursor and to the Brillouin precursor The 
rise of a (t) around t = Tg originates from the saddle 
points at 17^, which are then quasi degenerate and lo- 
cated in the vicinity of the pole at 17 = 0. The calcu- 
lation of the contribution ad{t) to a{t) of these three 
points requires to use an uniform asymptotic method 
[2J|. It is convenient to determine ad{t) through the cor- 
responding contribution hd{t) to the impulse response 
'^^^ exp [Q)] dn/2TT. Following the procedure 
we get hd{t) ~ be~'^*Ai[—b{t — Tg)] where 
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Ai {x) is the Airy function and b = (A^/St^ 
ad{t) = hd{x)dx reads 

ad{t) = e-''^s / e-'^^/''Ai {-x) dx 



,1/3 



Finally 



-T* / Ai(-a;)da; (9) 



the 2"^^ form holding when 7 <C & [2a|. ad(i) attains its 
absolute maximum at the first zero of Ai {—x) , that is 
for a; « 2.3 or i = i2 = Tg + 2.3/6. For h < t < t2, a{t) 
is well fitted by ai{t) + adit) (FigE]). For t > i2, is 
real and the frequencies are well separated (1^2 3> 7) . 
The contribution a2{t) of the two saddle points to a{t) 
can then again be derived from Eq.(l6]) with 62 = T^^/^^- 
It reads 



02(0 



-7* 



cos {n2t + AVgl^a/ i^l - A^) + 7r/4} 



fi2A 



(r!2 + A)- V (r!2 - A)- 



(10) 

The steepest descent contour T' passing through the four 
saddle points is now such that F + F' encircles the pole in 

= 0. The contributions ai{t) and a2{t) should then be 
completed by the corresponding residue, namely H{0) = 
e-''^". For t > t2, we get thus a{t) = e"''^' + ai{t) + 
02 (i). Again the agreement with the exact result is very 
good (FigdJ. The optical fields associated with 02 (i) 
may be considered as second precursors but, since they 
arrive after the rise of the main field , we suggest to name 
them postcursors. Contrary to those of the precursors, 
their instantaneous frequencies ujf (i) — ujc (t) are 
initially close to ujc before deviating from this frequency. 
Note that the oscillations in the falling tail of the pulses, 
observed in the experiments [iOl, are clearly related to 
our postcursors. 

We will now examine more briefiy the case of a medium 
with an electromagnetically induced transparency (FIT) 
window [13, m. 



29|. In such a medium, precursors have 
been indirectly demonstrated in an experiment of two- 
photon coincidence [s^l- We consider the simplest A ar- 
rangement with a resonant control field. If the coher- 
ence relaxation rate for the forbidden transition is small 
enough, the medium may be transparent at cue and its 
system function reads 



H (n) = cxp 



(11) 



where Vl^ is the modulus of the Rabi frequency of the 
coupling field (23,[2i|. We get then Tg = 2aL-f/nl. FigH 
shows the step responses a (t) obtained for different fir 
and for a value of aL intermediate between those of the 



celebrated experiments achieved on a lead vapor [27| and 
on an ultracold gas of atomic sodium As previously 
and for the same reasons, the very short term behavior 
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Figure 2: Same as Fig[T] for a medium with an electromag- 
netically induced transparency window. The parameters are 
aL = 600 and flrf-y = (a) 4.60 (b) 14.0 (c) 34.6, leading to (a) 
-yTg = 56.7 (b) -yTg = 6.12 and b = 1.397 (c) fTg = 1.00 and 
b — 4.647. Note that the group delays and thus the absolute 
time-scales are several orders larger than in the case of the 
natural frequency window. 



of a [t) (up to ti - 
In general, the 4*'* 



2qL7 



is given by Eq.©. 



degree equation giving the saddle 
point frequencies has no simple solutions but the follow- 
ing properties are easily demonstrated. Irrespective of 

0, fif [t) 



fir, S^2 i.'^g) — and, for t 



17 ± Q,r^jTg/ At 



while (t) ±inr/2. When < 7, Qt (0 and ^2 (t) 
keep non degenerate and purely imaginary at every time. 
If on the contrary fi^ > 7, these two frequencies coalesce 
at a time td > Tg in fid = i^r sin [sin~^ (j/Qr) /3] . For 
> 47, rid ~ il/i and td = Tg {l + A-y'^/inl) . ExpHcit 
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analytical expressions of a{t) can be obtained when ft,-/ J 
is moderate or large. 

In the first case, jTg = 2aLj'^/Q,l ^ 1 and the precur- 
sors will have a short duration compared to Tg. In this 
time domain {t) « i^t ±flr {1 + 3t/2Tg) ^Tgj^t and 



ai(t) 



TT 



7* 



cos [fJr (1 + V2Tg) ^ '^/^l 



1/2 



(12) 



If ^Tg is extremely large, the term t/2Tg may be ne- 
glected and ai(i) again equals a(^) given Eq.([3]). This 
particular case is examined in [3l|. When Sl^ < 7 or 
when fir > 7 with 7(irf — t^) » 1 (Figl2t), the only 
other significant contribution to a (t) is (t) associ- 
ated with the saddle point at fi^ (t) which tends to 
for t ^ Tg. We circumvent the difficulty due to the 
coincidence of the saddle point with a pole by pass- 
ing through the associated impulse response ft^ (t). It 



reads (<) 



with 9 



» -[nr{t-Tg)/4^ 

We finally get 



2n 

^ -^iTg/nl 
1 



(t) = -(l + crf [^Ir {t-Tg)/A^]) 



and 



where erf(x) is the error function. 



{t) 



(13) 



1 when 



(t-Tg) > 4 



and ai (t) + 03 (t) provides a good 



approximation of the exact step response at every time 

(FigEK). 

When fir ^ 7 the coupling field splits the origi- 
nal line in a doublet of lines approximately centered 
at iuc ± ilr/2. If, in addition, (17^/7)^ ^ 8aL/3 , 
then J {td — Tg) ^ 1 and the situation is analogous 
(but not identical) to that encountered with a natu- 
ral transparency window. The frequencies of the saddle 
points approximately equal f?^ « i7„ (t) ± fin (t) where 



In it) « (7/2) 



fin (t) is given by Eq. Q , with A 
ent contributions to a (t) then read 



(-l)"(l + 8Vr,)-^/' 



and where 
The differ- 



ai{t) 



a2{t) « 



rRe 



1 



[*(0+)— i7r/4] 



(14) 



Re<j _e[*K)+-/4] }. (15) 

\7r^" (fit) \ \nt ^ ^ ^ 



ad{t)~ I Ai (— .t) exp (— 7a;/36) (16) 



where b = {flj/12Tg 



,1/3 



As in the case of the natural 



frequency window, ai {t) + ad (t) and ai {t)+a2 (t) + H{0) 
fit very well the exact step response, respectively for 



< t < t2 = Tg+2.3/b and for t > t2 (Figl2jD and FigH:). 
The main difference is that a significant damping of the 
precursors is now compatible with a good transparency 
at ujQ- For intermediate values of fir it is so possible 
to observe both well developed precursors and postcur- 
sors without overlapping (FiglJb)- On the contrary, the 
tail of the precursors again partially interferes with the 
postcursors for very large rir(Figl2l:). 

To conclude, we have obtained, for the first time, fully 
analytic expressions of the entirety of the step response 
of linear media with a transparency window. Our re- 
sults show that these media, contrary to those generally 
considered, are well adapted to observe in a same exper- 
iment the precursors, the main field and the postcursors, 
all well distinguishable from each other and having com- 
parable amplitudes. Insofar as the parameters used in 
the calculations are representative of real experiments, 
we think that our work might stimulate an experimen- 
tal observation of these rich dynamics, which would, in 
turn, stimulate new theoretical investigations on related 
slow-light systems. 



[4 
[5 
[6 
[7 

[s: 

[9 
[10 

[11 
[12 

[13 

[14 
[15 
[16 

[17: 

[18 
[19 



Electronic address: 'bernard.segard® univ-lillel.fr 
A. Sommerfeld, Ann. Phys. (Leipzig) 44, 177 (1914). 
L. Brillouin, Ann. Phys. (Leipzig) 44^204 (1914). 
An English translation of [1J| and can be found in 
: L. Brillouin, Wave Propagation and Group Veloc- 
ity [AcAAemic Press, New York, 1960). 
J. A. Stratton, Electromagnetic Theory (McGraw-Hill, 
New York 1941) 

J.D. Jackson, Classical Electrodynamics, 2nd ed. (Wiley, 
New York 1975). 

K.E. Oughstun and G.C. Sherman, Pulse Propagation in 
Causal Dielectrics (Springer, Berlin 1994). 
K.E. Oughstun, Electromagnetic and Optical Pulse Prop- 
agation 1 (Springer, Berlin 2007), Ch.l. 
K.E. Oughstun, J. Opt. Soc. Am. A, 12, 1715 (1995). 
P. Pleshko and L Palocz, Phys. Rev. Lett. 22, 1201 
(1969). 

J. Aaviksoo, J. Kuhl, and K. Ploog, Phys. Rev. A 44, 
R5353 (1991). 

M.Sakai et ai, Phys.Rev.B 66, 033302 (2002). 

U. Osterberg, D. Andersson, and M. Lisak, Optics Com. 

277, 5 (2007) 

We use the definitions and sign conventions of 

A. Papoulis, The Fourier Integral and its applications (Mc 
Graw Hill, New York 1987) 

E. Varoquaux, G. A. Williams, and O. Avenel, Phys. Rev. 
B 34, 7617 (1986). See Eq.(49). 

B. Segard, J. Zemmouri, and B. Macke, Europhys. Lett. 
4, 47 (1987). See Fig.2. 

H. Jeong, A. M. C. Dawes, and D.J. Gauthier, Phys. Rev. 
Lett., 96, 143901 (2006). 

W.R. LeFew, S. Venakides, and D.J. Gauthier, e-print 
ArXiv 0705.4238v3. 

H. Tanaka et al, Phys. Rev. A 68, 053801 (2003). 
R.M. Camacho, M.V. Pack, and J.C. Howell, Phys. Rev. 



5 



A 73, 063812 (2006). 
[20] R.M. Camacho et al, Phys. Rev. Lett. 98, 153601 (2007). 
[21] M.R. Vanner, R.J. McLean, P. Hannaford, and A.M. 

Akulshin, J. Phys. B : At. MoL Opt. Phys. 41, 051004 

(2008). 

[22] B. Macke and B. Segard, Phys. Rev. A 73, 043802 (2006). 

[23] The results derived from Eq.(3) hold even if the two lines 
are not Lorentzian provided that they are so in their 
wings [20t]. The optical thickness aL/2 intervening in this 
equation may then be considerably larger than its actual 
value on resonance. 

[24] N. Bleistein and R.A. Handelsman, Asymptotic Expan- 
sions of Integrals (Dover, New York 1986). 

[25] V.A. Vasilev, M.Y. Kelbert, LA. Sazonov and LA. Cha- 



ban. Opt. Spectrosk. 64, 862 (1988) [Opt. Spectrosc. 64, 
513 (1988)]. 

[26] R.N. Shakhmuratov and J. Odeurs, Phys. Rev. A 77, 
033854 (2008). 

[27] A. Kasapi, M. Jain, G. Y. Yin, and S. E. Harris, Phys. 

Rev. Lett. 74, 2447 (1995). 
[28] L.V. Hau, S.E. Harris, Z. Dutton, and C.H. Behroozi, 

Nature 397, 594 (1999). 
[29] R.W. Boyd and D.J. Gauthier, in Progress in Optics, Vol. 

43 (Elsevier, Amsterdam 2002), p.497. 
[30] S. Du et al, Optics Letters 33, 2149 (2008). 
[31] H. Jeong and S. Du, Phys. Rev. A 79, 011802(R) (2009) 



